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II. — On the Theory of Planetary Disturbance. By the Rev. Bbice Bronwin. 



Bead November 30, 1848. 



1. JLN this paper I shall consider, with M. Hansen, the disturbance as af- 
fecting the radius vector (or rather the mean distance) and the mean longitude ; 
and I shall first, after his manner, employ two times. Various formulae not 
noticed by him (one of them fundamental) are given, in the hope that they 
may some time be made useful. The principal equations are investigated in 
a way that leads to some very elegant formulas, and the elimination of the 
quantities containing both the times is effected in a very simple manner. In 
finding the latitude, I propose to introduce the latitude itself and the reduction 
into the disturbance function ; by which means the part of that function de- 
pending on the inclination of the orbit to the fixed plane is greatly simplified, 
and the determination of the latitude and reduction rendered easier. But it is 
not my intention in the present paper to develope the functions in series of sines 
or cosines. 

The well-known differential equations of a planet's motion, referred to the 
plane of the orbit, are 



d 2 r A 2 u dR n 
df r» + r* + dr ' 



q dv , . (dR 
^dt= h=h -\~dv 



dL 



(1) 



p _ m'r m' _ xx' + yy' + zz 1 

K - -pt cos x- ( r '*_2rVcos x + <)* ; °° S X ~ rW ; 

where x, y, and z are the rectangular co-ordinates of the disturbed {z = 0), and 
x\ t/, and z 1 those of the disturbing body. The meaning of the other symbols 
is obvious. 
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When ml = 0, let the integrals of (1) be 

r o ~f (*> ^0' ^j *<>» *<>)> v o=<f> (*. «o» *o» *o« £ o)- 

And when the disturbing force is restored, let 

r =/ (<, a, e, w, e), w = (*, a, e, », e). 

Thus r and « are the same functions of t as when m' = 0, the elements a, e, &c., 

being variable, and determined in the usual manner. And, t being a new time, 

we make 

P =/ (t, a, e, v, e), \ = <p (t, a, e, ir, «). 

We further suppose, 

r t =f(z, a n e„ ir„ e,), v t =<f> {z, a„ e„ w„ e,) ; 

P,=f (£> a n e n *n 6 /)» \ = ( P(& °/» e n ■>• e /) J 
the elements a t , e n &c, and h n to be presently introduced, being constants, z a 
function of t, and £ a function of t and t. To these we must add the assumed 

relations, 

r = r,p, p = p^, v = v, + Gnt, \ — \ + Cn£. 

By changing t into f, we change p into r, \ into v, p, into r,, A, into v t , 
£ into /3, and £ into ■?. 

From what has been given above, we have necessarily rj —L=h l , 

P* -j- = h r And from the way in which a, e, &c, are found, 8r = 0, 8u = 0, 
the characteristic 8 denoting the variation relative to a, e, &c. ; therefore, 

But ^ = ^^' = ^tl = ^|- Consequently ^| = A ; or, 

3? "A?" () 



dt hp 2 h, 
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The two theorems now found are fundamental The latter is not no- 
ticed by M. Hansen ; the former, or one equivalent to it, is employed by 
him. But there is another fundamental theorem which he has not found ; for 

we have necessarily 

<P P tf « 

as in an undisturbed orbit ; because, since d operates on r only, the elements 
and h may vary in any manner whatever, and this equation will still hold. 

^- ^A%te_±_d h dZhd^d^_d^hd^d£ 

dr> ~ P 'di* "•" d$ dr dr h,? <# dr + A,£ d? dr ~ P < dr 2 + h£ d? dr' ° 7 ^ ' ' 

d*£ h 1 d* Pl , A 
= P'd? + ^d^ hytheBame ' 

But -j-j — j' + — = 0. Therefore, substituting for -73 and -7-^ tneu * values, 
ag p, p, dr 1 dg 

-_- — ~ and -,- — -4 , we have 
P* P* P, P, 



V-Hl-^, h * ( h l f>\. 

p> f-t'ds + h^Kft-py' 



or, since £p? = p 3 , 






p. (Pt; h* fl 

p % ~ Pl dr l Wg p) ; 



which easily reduces to 






There are still two other theorems to be investigated, which are given by 
M. Hansen ; but found by him in a manner very different from that which I 
shall employ. 

r = tf + F C0S (W ~ v) ^^ ~r ~f« = P C0S ( v ~ x ) = 

tie 

Tj (cos it cos t; + sin it sin «). 

e2 
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And differentiating relative to t, 

1 dr ue . , .dv . .dv, 

?dt = v Bm(v - 7r) dt> 0T > 6m< *^dt z=h ' 

X uf* UP lb(j 

T -T- = it sin (v — w) = y* (cos ir sin v — sin it cos u). 
A dt h 2 ' h 2 v 

Whence we easily deduce 

(1 u\ ldr . /xe 

^-pj cosu + A^ 8mt, = ^ C08,r - 

In like manner, 

t-p) c08X + ra sinX = F c08,r - 

/l fi\ . ldp fte . 

£. sin X - T -f cosX = £ sin -n. 

\p Ay hdr h 2 

U& UiB 

Equating the two sets of values of p cos ir and p sin it, 

(1 fi\ ldr . /I A ^ ldp . . 

{r-£>) c ™ v + kdt amV = {}-h>r aX + hdr SmX > 
(\ M \ . ldr (1 /A . . ldp . 

{r-I 2 ) SmV -hdt COaV = {p-h 2 ) SmX -hTr C09X - 

Multiplying these by cos v, sin v, cos X, and sin X, and adding and sub- 
tracting results, 

\-¥ = (cp-^ C °^ X - V) + lfr S]n ^- V ^ 

ldr /l it\ . „ .1 dp .. . 

hdt=-{-p-h 2 ) sm ^ x - v ^hi C0 ^ x - v ^ i 

1 fi (1 /A . 1 dr . 

p-F = V7-I 2 J co8(x - w >-A^ 8in ( x - t;) ' 

1 dp /l ,A . /x x 1 dr 

ra = (r-Fj 8m ( X -^ + A^ C0S ( X - V )- 



(5) 
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The equations (5) may be put under the following form : 

-j (r sin v) — -j- (p sin X) = r (cos » — cos \). 
dt At n 

57 (r cos w) — -r- (p cos \) = ~ (sin A — sin v). 

dt WT fl 

The first and second of (6) may be written thus : 

-=^|l-cos(\-«)}+^^|/>sin(\-«)S; 
dr p . ,. i . d , ... ., 



dt h d 



T 



These, and the third and fourth of (6), are very elegant equations, and are 
deserving of notice on the ground that, possibly, some use may be made of 
them. I shall only employ the third of (6). Differentiating it relative to t, 
there results, 

1 dp C2ti 2« /A .Idh fld\ ,_ N 1 dh . ,. . 1 dr 

7l = {i?-Y cos ^-^h + {hdt co8 ^-^-¥dt 8m ^- v n-dt 

Ehminating -7- by the third and fourth of (6), and substituting for -7-j — 3 

+ 4 lts value from (1), we find 

I dp IdpdX fl u /l u\ ,„ N 1dA 1 . /x . rfi? 
?i = Ar T ^ + b + ^-(r + FJ C0S ( X - u) ]^-A Sm(X - w) ^- 

Multiply this by h, put for -7- its value from (1), and for -^ in the first 

member its value -5- ; then making 



we have 



MH)^-»>-(H)}§-°^-<. 
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dp d\ dp d\ 



dt dr dr dt 

which is a very remarkable equation. 
From this we deduce successively, 



= T, (7) 



dp d\ _dp d\ _ T ^ dp 
dt dr dr dt dr' 

K(±d£_dpd£\_ dp 

p?\dtdr dr dt) - + dr' 
but 

dt~ Pt dC*d£ dt 1 dr p 'dr* d£ dr' 



Therefore, by substitution in the above, 

h 
P. 



This again is a very remarkable equation, and has a remarkable corres- 
pondence with (7). Dividing it by — -^ = —3 = — , we have 

/>/<*»■ Pj? pt 

l d I\ 

d l^ d l * _I ,7.,^ d ± 
dt dr\d£\ ~h^ + h P dr'' 

or, 

d * - x ^rJ nk d W j. d * [51 n\ 

Jt~h ( ^ I+ 'U~dr~ + dr\d^ ( ^ 

\Trl 

Equation (7) multiplied by 2p gives 

But differentiating p 2 -5- = A, relative to *, we have 

, d*\ dpd\_dh 

f drdt + P dtdr~lt' 



or, 
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The above, subtracted from this, gives 

This is singular, and very worthy of notice. This theory, as we have 
treated it, gives rise to many very interesting formulae. 

Therefore (10) becomes 

\(K* d s\±t„m^ 

dt 



Or, since A^ 2 =^-, 



dr \ /S dt) dr "* ' 




2. In the preceding section we have found all the fundamental formulas, 
and indeed more than are absolutely required. We now proceed to discuss 
some points preparatory to integration. 

We might immediately integrate (11) relative to t; but the integral of the 
second member would contain a very great number of terms with r in their 
coefficients, or of the form rf(t). We must, therefore, proceed otherwise. 

Make g = t + m, £ = t + 0. Then $ = J^, d^ = 1+ ~d' Substitutin g 

these values in (9) and (11), neglecting terms involving the fourth power of 
the disturbing force, and putting 



h r 2A dr h dt h r h dr 



they become 
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dji_ p d£ dj> _d§ # d$ 

dt ~ dr dt dr dt dr' 

^.-9 ^ (dj> dj> _dj> dtp? 

dnti~ dr\dt dr dt dr* 

Let the integrals of these relative to t be 

<l<p 



) 



(1) 



* =/(t) + *(t,0, f T =/,(t) + «(t,0 : 

J\r) and /,(t) being the arbitraries of the integration. If we make m' = ; 
then g = 0, and \ = \ . But in this case X is a function of r without t Con- 
sequently, \„ £, and 0, are functions of t without * ; and, therefore, -y- = 0. 
Hence, all the terms in the second members of (1) vanish, and we have 

| =/( T ) ( ^ = / /(T ); 

£ and -~ being the values of £ and -j- when m' = 0. 
«t ° dr 

Now if we make the disturbing force to vanish in (2) and (4) of the first 

section, they become 

d£o __ h 
dr 



htf 
d%_]' L (h 2 

1+ dr ~htf' 



(2) 



the last of these being derived from the first by putting for £ its value. These 

equations, being integrated, give us £ and £ , or & and fa. 

h 2 
To integrate the second of these, make £ = ~ (1 +7), when it becomes 

,,3 

To abridge, make ^- 3 = l+eE 1 +e 2 E 2 + &c, where E„ E u &c, are known 
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functions of the cosine of w»t + Co — t and its multiples ; and assume 

y = B + eB 1 + e*B i + .... 

By substitution, and equalling separately to nothing the terms multiplied by 
the different powers of e, we have 

d*B .„ _ (PBi .„ „ <PB . . 
-^ + nJB = 0, -^5! + n 2 B 1 +E l _ = 0, &c. 

The first of these gives B = c cos (w„t + k). This value substituted in the 
second, it will give Pi ; and so on. But in reducing the periodicals to one ar- 
gument, we should have a series of sines and another of cosines, unless we make 
k = to — ir . In this case we must have % — w„ = e y — n /t and making n = 
n, (1 + 6), we must develope thus : 

cos (n„T + en — 7r ) = cos (n/r + e„ — w ) — bn/r sin (n y T + % — v ) — &c. 

The use of this is, by suitably determining b, to take away an improper 
term. In some cases, perhaps, both sines and cosines may be necessary. 
But perhaps this will be more conveniently done thus : 

- = p {1 + <?o cos (\> -<*,)}, — =-£• |1 + <?, cos (X„o -*,)}. 

But \, j0 = \. Therefore we must have x = *•„ or we should have both 
sines and cosines. And then we shall have 

( h -l - lV = cos (\ - *•„) = cos (\„ - W/ ) = fA. _ i)l . 

or, after a little reduction, 

^1 — 1 _ ^ ^ ^ 
W»~ e , W»o e ,' 

But po = p^oZo. Substituting this value, we easily find 

or, x " ' 

l_ e / -g» />„ a/, (!-«») 
& 1 - « 2 Oo^/ <V/(1 —4?) ' 

VOL. XXII. p 
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Make «b = e, (1 — a), a being a small quantity of the order of the disturbing 
force ; and we may write 

So <h Oo 

To which we may add from the third of (2) ; 

d<po _ Oo* C_ _ .. 
dT~a~*tf 

A, B, and C being functions of e 2 , and a easily found. But /» /j0 , a function of 
\, and, therefore, a function of tv, must be converted into a function of n,r 
by Taylor's theorem. The results obtained are much more simple than those 
obtained by M. Hansen in the same case, and I think also more convenient. 
If we wish to make * = *, + V> we have 

cos (\, — * ) = cos 17 cos (\ — 7r ) + sin rj sin (\, — » ). 
— -p = sin (\ — •*<>) = cos r\ sin (\ — ir t ) — sun/ cos (X — ir t ), 



cos 



(\,o ~ *,) = cos (\ - *,) = (A- - 1 ) \ . 



Between these we may eliminate sin (\ — «•,) and cos (\ — * /} ) and obtain 

a result involving -p ; but I shall not pursue the subject further. 

We may satisfy (2 ), and all the requisite conditions, by simply making e = e„ 

» c = "■/» £o = ~ ? » an ^ w ,£o= V J but tms would not leave us any arbitrary con- 
es, 

stant except a , and we might have an unsuitable term which we could not 

get rid of. 

We must now proceed to another class of formulas, some of which will be 

wanted. But I shall not confine myself to these. 

Let Vj-U ( jt-J denote that r is to be changed into t in these quantities, 

after the operation of differentiation is performed. And thus this change will 
be denoted in other cases. Thus from (2) of the first section we have 
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ms_± 
(d±\_ _h__ 

\drj-kp l - 



(3) 



Eliminating r— 5 between the first of these and (3) of the first section, there 

results 



\dr) ~dt + h, ' 

or, 

fd<j>\ _ dw Gnr* ( 

\d^)~di + 17" 

But it is obvious that f-^j + f-^j = -,- . Therefore, 



(4) 



6nr* 

IT' 



(5) 



But f ^7 J =^= 0, from the manner in which the variable elements are 
found ; and -^ = -^-. Therefore, 

which by (5) gives 

and since also we necessarily have 

\drj + \dtj ~ dt ' 
F2 
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therefore, 

\dr)~dt \ dz' l'> 

The preceding are all the formulas of this kind that are really wanted ; but 
I shall give some others, on the ground that they may possibly be some time 
found useful. If we could change t into t in higher differentials than the first, 
we might, perhaps, by differentiating the fundamental theorems, or any of the 
others, relative to t or t, and then changing t into t, find some new theorems 
which might be made available for simplification, or some way useful. 

From (2), section (1), or -3- =-rj 2 - 1, putting r + for £, we have 
djl__2hdj Therefore 

dT>~ hjedr' ineretore ' 



or by (7), 



<fr* j ~ hjf \dr) ; 



V***/ Kjp dt+htf dz ' K } 

Equation (4), section (1), gives immediately 

And from (2), referred to above, after putting r+ tj> for £, and differentiating 
relative to t, we have 

££_J_ dh_ 2h dt 

drdt~h l t?dt kg dV 



/<*>\ 1 dh 2h (d£ 

[drdt) ~ h,p dt hfi \dt 

or, by substitution from (6), 

1 dh Snh d(r*) 



3^\ 

\drdtj 



hjp dt AJ/3P dz ' ( 10 ) 

I might extend this list, but not without some difficulty ; and as I am not 
sure of its utility, I shall here leave it. 
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We now turn to the equations (1) of this section. The first contains both 
£ and <f> ; the second only <p. This last, therefore, is preferable ; and as we 
do not want them both, we shall take this, which reduces to 

drdt ~° + drdt Wt dr 1 r dt dr* V dr 



For the first and second power of the disturbing force only, 

d'<j> _ „ d<(> dty d<f> d*<f> _ „ Q d(p d<f> d? <f> 
drdt ~ dr drdt dt dr* ~ dr dt dr 1 ' 

Putting this value in the second member of the above, we have 

d^<f> _ „ _ d<f> d<f> <^0 dtp d<f> <P<p 
drdt dr dt dr* dt dr dr* ' , 

c r 2fi P /2 2/A _ ."I IdR . ,. .2 P dR 6nd( P 2 )t ( ' 

3. The values of <p and £ would be very troublesome to find ; for they 
would contain a great number of terms having r — t in their coefficients, and 
which in w and /3 would vanish ; and many which, when t is changed into t, 
would unite with others. If, therefore, we could get rid of these quantities, 
we should greatly diminish the labour of integration. This, happily, we are 
able to effect. To accomplish it we will develope S relative to these quantities. 

Equation (2) of section (1) gives S = » (l + ^= * (l - *g + 1*£). 

Also d( ^ - rf( ^> - * rf W> 4- o^^ - * d( ^ * 4- „> d( ^ k d W> a. 
A1S °' ~dr dT-S ~dr~ + P '~dr~ ~ ^ ~dYTr + P '~dT = h~af + 

p* \ — . These values are to be substituted in that of S, after we have 
dr 

changed p into p£, v into v, + int, and \ into \ + int. Make, therefore, 
/ox fAJ.2/v, /2 2/A ."I 1 dR , . ,_ , 2/»,rfi? 

^i£> and 

A. dr * an(1 
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,_ f2/u/>, /2 2fi\ ,. .1 1 <£R . ,. .2/., di? 

^={^-{r + i)^ cos ^'^)m^ +sm(x '- v ' ) 4^ = 

*°' h dv + h, d£ ' 
Putting the above values of £ and -A— in that of &, after developing rela- 
tive to £ ; we have, making /»* the value of p, (a function of t + 0) when is 
made nothing. 

S-(S)-(T) (£_ -^j -_,>,_ _ 

We have now got rid of £. Let (S°) and (T ) be the values of (S) and 
( T) when is made nothing, or £ = t ; then, by Taylor's theorem, 

(S) = (S .) + ^ + W>, (2) 

neglecting higher powers of the disturbing force ; since we exclude those of the 
fourth order. Now <f> enters nowhere but where it appears. 

And putting — for - -v- for convenience, 

d(T°) = d(S°) en d?(p™) 

dr dr hj di* 

Therefore 

Put in (1) the value of (S) given by (2), and that of (T 7 ) given by (3) ; 
and we have 
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'-o(w^i3)**«P('-»$-iS*(i-iy 

dr dr \ * dr) A, dr 2 \ drj A, dr dr 2 



2A, 



2 dr 8 2A, cM <2t ' 



and 



° rfx ~ ^ ' V^t * *»/ + 9 dr dr 2A dt dr 2 

Cn_ d( P ?) dtf Gntf <fy dfy 
2A, dr dr 2 + A, dr dr 2 ' 

We must now put these values of «S and S-j- in (11) of the last section, 
and we shall find as the result, 

d+dtyt \_df\_in d( P f)d<pf d<t>\ tnpj o?*/ _ d<j>\ <f> 2 d\S°) 

dt dr 2 V dr) 2A, dr dr \ + * dr,/ + A, dr 2 V dr) + 2 tfr 2 

en d 2 ( P ?) <pd<p Smp d(ff) d*<p 

A, dr 2 2 dr h t dr dr 2 ' K } 

We can do nothing with the last in its present form, on account of the 

terms which contain <f> and ■— . But happily these may be driven out by a little 

transformation. Neglecting terms containing the third power of the disturbing 
force, (4) gives 

(cox _ d3 <P K&\ d< P A d ( S °) d< t> <*V 1 dh d<f> en P <>?d 2 <p 
V ' ~ drdt * ^ > dr * dr dr dr 2 + 2A dt dt ~ h, 1? 



6n d( P 02 ) dtp _ 
+ 2A, dr dr ~ 



d 2 <p 

drdt 



* dr drdt * dr'dt dt dr 2 + 2A dt dr A, dr 2 + 2A, dr dr' 



With this value we eliminate («S°) from (4), in those terms where the dis- 
turbing force rises above the first power, and we thus find 
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This has been put under the simplest form. It still contains three terms 
which we cannot manage ; but these may be easily eliminated. Neglecting the 
powers of the disturbing force above the first, and integrating relative to t, add- 
ing the correction -^2, or the value of ^— when the disturbing force is made 

to vanish, we have 

& = ^» + J (S°) dt = W suppose. ( 6 ) 

Now if we include both the first and second powers of the disturbing force, 
we have 

which gives us immediately, 

d<f> 



dr 
And 



d<p 2 .d'<p_ w rd<f>dh Gnpfrdty 6w d(p°?) [dtp 

d?~ <p d^- yr ~^}TrT + ~X\'aV dt ~2h t dr Idv'- 

't 2 ^ dr dr* ^dr^dAdr i dr 2 dr> ]' 






Make 



d<p .dp d 2 (p Y 



Now as ail the terms in the value of X after W are of the second order, 
and will be multiplied in (5) by quantities of the first order, we may put W 

for -j- in the value of this quantity. Then 
and (5) will become 
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a^di- {b) + dt\ i dT> +(l> dr> *dS VVdrdi* *+ <h») 
1 dh Gn dip?) Y faff dX 
2h*dt~2h, dr ^ + h, dr' w 

J ■ 2 J2 J. 

We must now eliminate -^ , -y?, &c. For this purpose we take 
which squared, neglecting terms above the third order, gives 

dW 

dr 

dW 

- W* + 2WX+ \W> + 2<pW^-. 
dty _dX d<p dty ,(P<p_ 

di*' dT + *dr dT* +<f> d?~ 

dX , «„ r dW &W 

Substituting these values in (8), it becomes 

— x—-— d il¥)x en ff ^ 

2A ^ 2A, dr + A, rf T ' 
Let 

Y=-lW>+lWX-VX^-^ d ^ i Xdt + *f\ d £dt. (9) 
By means of this the integral of the preceding is 

£-,r + ,(£ +i <*) + **£ + r. do 

VOL. XXII. Q 



42 The Rev. Bkick Bbonwtn on the Theory of Planetary Disturbance. 



Since ( -p ) = -=- + —r 1 ; changing t into t, we have 



dw 
~dt 



= ( ^)__, +M ,^_ + ^^ r j + K(-^)+(n (11) 

We have now got entirely quit of £ and ; and t is contained only in p° t and 
\°, which are the same functions of r and a n e n &c., as in an undisturbed orbit ; 
and we may everywhere put them without the sign of integration, and then 
change t into t. Thus we have virtually got rid of t, and the result may be 
put under the same form as if it were obtained with only the ordinary time t. 

If we change z into t + w, (3) of section (1) gives 

h _ 1 dw Gnr*, 

hj*~ + ~di + ~i;' 

and 



a h* f , . (dw 6m*\ „ "| 



which will give /3 when -j- is known. And we may change /3 into e*, and take 
the logarithm of both members, if we wish to have a result in the experimental 

a 

form, which is that of M. Hansen. Or we may change /3 into 1 + — ; then 

we have 

r = r, + /3. 

But, changing £ into r + (f>, we may employ (2) of section (1), which gives 
and 






If we wish to have a result in the experimental form, this is, perhaps, the 
most convenient. And thus changing £ into e*, and taking the logarithm of 
both members, we have 



= iO-£>og(i + £). 
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M. Hansen employs this method, and differentiates relative to t, by which 

means he gets rid of £ log ( t- ) • Thus we have 






di* 



1 + 



d<f>' 



and changing r into t, 



3)~* 



^0 



,1 + 



rf0 



Or, by (7) of the section referred to, 



dfi _ tnr dr, . 

dt ~~ k, dz * 



dfy 
d^ 



°L± \ 



1 + 



rf0 



In this form M. Hansen leaves it. But we might employ (8) of the same 
section to reduce it, or to give a distinct form ; but that would introduce 

h 

r- again. 



4. We now proceed to the determination of the latitude, and the reduction 
to a fixed plane. Make a the sine of the latitude, i the inclination, and S and 
the longitude of the node on the plane of the orbit and on the fixed plane 
respectively, $ having an origin fixed on the former plane. By making the 
plane of the orbit to turn round the radius vector an infinitesimal space, it is 
easily seen that, 

d$ = cos ide ; (1) 

and we have the known formulae, 

di _ 1 dR _ cos i dB 



h sin i d0 h sin i d§ ' 
d$ _ cos i dR 
dt h sin i di 
g2 



(2) 
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I shall give some theorems here similar to (6) of the first section, without, 
however, making any use of them. 

We have o = sin i sin (v — $) = sin i (sin v cos 3 — cos v sin 3), -7- = sin i 

ax 

cos (v — 3) -5- = — 2 sin i eos (v — 3). Therefore, -r -r — sin i cos (t; — 3) = 

sin i (cos v cos 3 + sin v sin 3). From these we easily find 

1* do . . 

a sin v + -7- — cos = sin z cos 3. 
A at 

r 2 do . . . . . 

— <»• cos v + 7- t- sin w = sin z sin 3-. 
ft at 

Again, let «• be the same function of t and the variable elements that o is 
of t and the same elements ; so that 



k = sin i sin (A — 3-), -j- = -^ sin i cos (\ — 3); 



C?T p' 

and, as before, we shall have 

. . p 2 die . . . . 
*• sin \ + V -7- cos \ = sin 1 cos j, 
A dt 

p 2 die . . . . _ 

— «■ cos a + V -7- sin X = sin * sin 3. 

Equating the two sets of values of sin i cos 3 and sin i sin 3-, 

r 2 da . _ o 2 die 

a sin w + -7- -r cos w = *■ sm a + -=- -7- cos A, 
A dt n ar 

r 2 do . - p 2 die . 

o cos v — - r - r 8in« = «f cos a — r -7- sm A. 
h dt h. or 

By multiplying these by sin v, cos t>, sin \, cos \, and adding and subtracting 
the products, we shall have no difficulty in deducing 

= k cos (\ — v) ^- — -7- sin (X — u). 

r 2 rfff . ,_ . p 2 die , x . 

I3 =««n(X-.)+J s eoi(\-.> 
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,, . r 2 da . ,. . 
K = a cos (A — v) + -r -3- sin (A — V), 

p 2 a"* . /A N , r 2 rfff / . N 

^ -r- = — <r sm (A — v) + -7- t, cos (\ — w). 
h dr v h dt 

These two sets of equations correspond to (5) and (6) of the first section. 

_ r 2 do . . . . . c?<7 . . . -v , rfff 

From j T = sin« cos (u - 3), -jt = - sm i cos (v — 3), and -p = cos 1 

sin (w — 3) ; we easily find 



<f<r r 2 da da cos i 



= a 



d$ h dt ' rfi sin i ' 

die p 1 die die cos i 



= Af 



d3 h dr' di sin i 



(3) 



which may be found useful. 

The part of the disturbance function depending on the inclination of the 
orbit to the fixed plane is very troublesome to express by means of i, 3, and 6, 
and the corresponding quantities i\ 3', and &, relative to the disturbing body ; 
and when either 3- is expressed by means of 0, 3' by means of (f, or the latter 
by means of the former, would contain a great number of terms. But it may 
be expressed very simply without these quantities, by means of the latitude <f> 
and the reduction A, and the corresponding quantities <j> and A' relative to 
the disturbing body. Thus we should have v — A, and xf — A', and v' — A' for 
the longitudes on the fixed plane, and 

tJC 7/ 21 

- = cos <f> cos (v — A), £ = cos <f> sin (v — A), - = sin <f> , 

stf y f g r 

— = cos eft cos (1/ — A'), —, = cos 0' sin (»' — A'), — = sin ep'. 

„ . xx' + yy' + zz' . „ 

But cos x = / ; therefore, 

rr 

cos x = cos eft cos 4/ {cos (v- A) cos (v'~ A') + sin (v- A) sin (v' - A') J + 

sin sin 0' = cos cos ef>' cos (u - */ - A + A') + sin sin 0' = (1- <r 2 )* 

(l-<x' 2 )» j cos(v_v')+sin(v-v')(A- A')-cos(»--»')i^=^- - &c ]+«/. (4) 
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We suppose here that a and </, and therefore also A and A', are small, as in 
the case of the moon and the old planets. 

Making for a moment r t the projection of the radius vector, and v t the 
longitude on the fixed plane, and comparing the differentials of the areas des- 
cribed on this plane and that of the orbit, we have 

cos i r*dv = r* <fo ; = r* cos 2 <f> dv n 

and cos i dv = cos 1 <f> dv t = cos 2 $ (dv — dA) ; or 

(cos z 0— cos2)-j7=cos z 0-^-; or ( 2 sin 2 ~ — sm 2 0] -j=co8 2 0-p. 

Whence we derive 

o • 2 * • 2 ^ . r2 ^ 2 ^ 
2 sm 2 - = sm 2 4>+-r -jj cos 2 <f>, 

dA h U • 2 »' • 2 J\ 

And putting a 2 for sin 2 0, 1 — a 2 for cos 2 0, 



c ■,» , , r 2 dA „ 

rfA A A> • 2 » j\ 

^ = ?"(l37)( v 2sm 2-^- 



The equation iS=cos id0 gives rf0= . dS, and d (0 — 3) = 

COS % 

2sin 2 ^ 
</ (0 — 3) = — c&. And by integration 

COS % 



COS I 



(5) 



1 ) d$ ; or, 



f sm 2 



— 3- = 2 j r dS = a suppose. 



COS l 



(6) 



Let A, = (u - 3-) - (v, - 6), where v, = v - A. Then 

sin A, = sin (u — 3-) cos (y i — 6) — cos (v — 3) sin (v / — 0). 
But by the well-known theorems of spherical trigonometry, sin (v — S) = 
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s ™4 t sin ( . _ e) = *^t, and cos («,, - 6) = 22?lizl) . Therefore, by sub- 
sin i K ' ' tan i \ i / cos ^ 

stitution, 

A — tan ih oe\a (m 5^ I „„,„ 

^sin i tan ij 2 

sin i . , Q x 

■. tan cos (v — 3). 



sin A, = tan A cos (» - 3) ( - — ; . J = tan - tan <b cos (v — 3) = 

* r \ \ 01T1 _ ton it V ' x 



2 cos 2 J 

2 7 

But sin A, = A, — %A*, and sin i cos (v — 3-) = j -y-. Consequently, ,A — 
3 tan r 2 <2<r _ r 2 <r da __ r 2 rf \/(l — it 2 ) 

6 Z 2 

A a r 2 dS(l-a*) l1a3 

2rtcos S g 
Now A, = » — » y — 3 + = A — 3- + = A + a ; and, therefore, 

2 A cos 2 ^ 

neglecting smaller quantities, observing that a is of the order of the disturbing 
force multiplied by sin 2 - . 

We want now to find -j- and -7^. To do this, we shall transpose £A 3 , 

and, after taking the partial differential, divide by 1 — J A 2 , or multiply by 

, , 1 .. j r n d 2 a h , cos i da , d 2 a 

1 + £ A 2 ; and by means of -p-j- = -3 cos « cos (v - 3) = -, — : -3- , and y— r = 

didt r* v ' sin % dt d$dt 

-3 a, we shall find 

dA r 2 da a f . /2 - a i \ il , , , v 
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9 ' 2 ' 

^ o ,«/(l-'r)\ h i d#l-o i J K * ' cos e v * 

a COS ^r 

r* d*? 
But p -p- = sin 2 i cos 2 (y — $) = sin 2 i — sin 2 i sin 2 (» — 3) = sin 2 i — a 2 . 

Therefore, putting this value in the preceding, we have 

But the formulae we have obtained are not convenient for actual applica- 
tion ; it may be well, therefore, to give them in series. For this purpose we 
make sin i = s, sin (v — 3) = rj ; then a = sr/. But as, in the first differentials, 

the elements do not vary, we must make -j- = s -r, and in — we must not dif- 

J dt dV dt 

ferentiate §. 

Now, expanding */{\ — a 2 ) in (7), we find 

2Acos 2 - 
z 

neglecting smaller terms. Or, 

2Acos 2 5 
But 

1 2 2 1 — cos i 1— v(l — s 3 ) . , , , , 
._ — _ _ \_ L-l + ^f + Jgti 1 . 

a . l . . 4 % . z s 2 s 2 s 2 ^ a io 

2 cos 2 - 4 sin 2 - cos 2 - 
2 2 2 

Substituting this value, we have at length, 
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where some small quantities are retained, which, perhaps, are not wanted. 

In taking the partial differential coefficients of A, it will be better to find 

dA , dA 
-j-, and not -jv. 
ds di 

Write this for a moment, -r- = M + ^A* -^-. Transpose the last term of the 

second member. Then (1 - £a 2 ) -t- = M; whence -r- = M (1 - |A*) _1 = 
^(1 + ^A 2 ); and 

^ = TI^" + ^'' + ^ +5, ^ + ^ 5 + ^Kl + ^). (11) 

Preparatory to finding -j- , we may observe that -^ = cos (v — $) -^ = — 

cos (v - a) ; and, therefore, -^ = -j sin (u - a) = -5 »;. Also, -£ = - cos 

(w-S) = -^; and^^J = cos 2 (tt-$) = l-sin 2 (w-$) = l-^. We 
now find 

dA 2sia2 i 

'(i-^ll + ^rt + W + ^dV + AV+MV)}. 

By further reduction, and treating the equation with regard to A as in the 

last case, 

2 sin*i 

( 1 + * A2 >-^sT( 1 + ^>- 
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But l!!l! = l^2ii= _L_ 1= *__ 1 = ^ + fs <_^« 

cos« cos i cose v (1 — «■*) z a lb 

Therefore, 

<ZA 



^=^-i+^+« 2 (-|+K+'/ 4 )+^(-f-TW-W+K)}(i + 4^). (12) 



Since 2 sin 2 ^ = 1 - cos i = 1 - /(l - s 2 ) = £s* + %s* + fa*. Putting this 
value in (5), and developing the terms, it becomes 

^^^(h-^+^a + W-^ + s'^+W-rhrf-rf)]. (13) 
We will now transform (2), so as to introduce a and A, or rather rj and A. 

ds _ . di _ cos 2 i <£R _ 1 — s 2 dR 

dt dt h sin i d$ hs d$ ' 

^ - __ £21*2 ^? - _ 1 ~ s2 ^ 
cfr A sin i ds ~~ As <fo ' 

Or, putting s»; for <r in 5, and introducing the partial differentials of the new 
quantities, 

ds _ 1-g 2 (dR dn dR dA\ ] 

dt~ hs \drj d$ + dA d$)' | 

d$ _ 1-s 2 (dR dRdA\ i ^ 

<fo ~ hs \ ds dA ds }' J 

It may not be amiss to find 2 sin 2 - separately. Thus, i being the mean 

value of i, or its value when the disturbing force is nothing, and, therefore, a 
constant quantity; we have 

2 sin 2 J = 2sin 2 § + 2 Jdsin 2 £=2 sin 2 J + J sin idi = 
. dR _ a ._ 2 i _ ( dt . (dR dt) dR dA 



_ . J , (dt .dR a . Jo , (dt 



dvi dS dA d$ 
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„ . 1 dR dfj 1 dR dA n . , . , „,, 
Make O,-J + AdA^ =P ' t0abndge ' Then 

2 sin 2 | =2 sin 2 5° + Jcos iPdt = 2 sin 2 £ + \Pdt - 2 Jsin 2 £P<&. 
2 2 A & 

But substituting this value of 2 sin" ^ under the integral sign of the second 
member, 

2 sin 2 % - = 2 sin 2 \ + cos 4 J i^* - £ (f-P*) 2 + 2 !•?* I sin2 5 m 

By continuing these operations, we find 
2sin 2 | = 2sin 2 | + cosi,jP^-^cosio(jP^) 2 + ^3 COsi o(J^0 3 - &c - ( 15 ) 

This may serve to eliminate sin 2 -. But since 2 sin 2 - = 1 — cos i = 1 — 

V^l — sin 2 i) = £ sin 2 i + £ sin 4 i + -jJ ff sin 6 i + &c., it may also serve to elimi- 
nate sin 2 i and sin i, if we should find it convenient to do so. 

It is much easier to find r and v on the plane of the orbit, than to find the 
values of the corresponding quantities on the fixed plane ; but it is more diffi- 
cult to find the latitude in the former case. It is very troublesome to find it 
directly from the variable values of the elements i and 3- ; and yet we must 
have the values of these quantities separately to a considerable degree of ex- 
actness in order to find a. M. Hansen has found the latitude by means of sin i 
sin 9, sin i cos S ; which is a much better method ; and he has found a by means 
of these latter quantities, or rather functions derived from them. But this is 
attended with a great deal of trouble. I propose to pursue a different course, 
and to find sin i, or s, separately, by the equation given above for the purpose , 
and then to find rj. To do which I shall find 

y = sin (x — 3 — <ynt), 

x being a constant quantity. This is finding a function of 9 + ynt. 

dt d$\^^dt) 

dy f I-* 2 (dR dRdA\\ dy „ 

d^{^--hT-{-dJ + dA-dl)} = I Qm ^ SQ - 

h2 
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It must be observed, that ynt is the uniform regression of the node, and y is to 
be so determined as to take away the constant terms from Q. 

Now y = sin (x — %) the value when the disturbing force is made nothing . 
therefore, integrating, 

y = sin (x - 3 ) + J^ Qdt. 

For a first approximation, we make -^ = — cos (x — 3 ) ; and after the in- 
tegration is performed, we may make x anything we please. We shall make 
x = v + ynt = w / + tnt + ynt. Then y = a, and we shall have 

a = sin (0 + ynt - 3 ) + f ^| Qdt. (16) 

After this substitution has been made for x, and v replaced by v,+ Gnt, v Jt being 
a given function of z, may be allowed to remain, or may be developed in terms 
of z, and z may be developed in terms of t. 

"We may make y = sin i sin (x — 3 — ynt) = sin i cos (3 + <ynf ) sin x — sin 1 
sin (3 + yraJ) cos x =p sin x — q cos ar. 

p = sin » cos (3- + ynt), q = sin i sin (3- + ynt). 

dp . .di . . . , n . (dSr \ 

-j- = cosi cos (3 + ynt) j — sin 1 sin (3 + ynt) I -7- + yn J . 

* = cos i sin (3 + ynf) -7 • + sin i cos (3 + ynt) (-7- + yn). 

Or, 

J» cos i di fd& \ p ds /rf3 , 

tto cos i <ft /rf3 \ 9 <fo , /d3 

i = ^sini^ + ^U +7n y = ^ +;> U +7W 

We may substitute in this, for — and -7-, their values from (14) ; but as I pre- 
fer the former method, I shall not pursue this further ; nor is it necessary, since 
any one who is desirous of doing it may easily carry it through. We may ob- 
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serve, however, that making Ip, Zq, 8s, 83-, the alterations produced in p, q, &c, 
by the disturbing force, we have for the first power of that force, 

hp=£ss-qi$, hq = ^is+pl^; 
s s 

whence, observing that s 3 = p 2 + q 2 , we have 

s ' s 2 

Thus we find the alterations produced in s and 3- from those produced in p 
and q. And in like manner we may find the alterations depending on the se- 
cond power of the disturbing force. 

In this section we have taken no notice of the development of r in the for- 
mulae where it has appeared. Making r=r t /3, we may either let it stand thus, 

or put from (3) section (1) its value in terms of -y-, which we easily find to 
be 






Then we shall have only r, and v,, besides the terms depending on the latitude, 
to develope. These are given functions of z = t + wj, and may, therefore, be de- 
veloped together by Taylor's theorem. 

We might have developed the value of S in (11), section (2), relative to £ 
and 0, differently. 

Since the difference only of \ and v enters into the composition of S, we may 
change these quantities into \ and v t , or into the anomalies \ 7 — •n l and «, — »,. 
Using for simplicity A, and v t , and making 

,,_ (\ fi\ 2 cos v t dR 2 sin v, dE 
\r h 2 J h dv h dr" 1 

_ /l fi\ 2 sin v, dR 2 cos v, dR 

~ \r + h 2 ) h ~fo h W' 

we have 
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„ A 2fi P \ ldB % ._ . . , r e» d( P 3 ) 

s = { 1 + -W)hd;-t ,cosX > M -<' amX < N -T dT 

We eliminate | from this by putting p, -j-\ ( 1 — \ -? + | ~o ) for p = />,£. 

A* 
Let (>S) be the value of aS when p is changed into p, t-j i and we have 

Gn ^ r 2 /<fy_#^\l 
A, <*rl/'V*r *»y]' 

Eliminate ( P/ cos \M + Pl sin \JV) — U?- § -5. j by 
and we find as the result, 

Make (P) = (5) - J ^ + £ S£l. Then, 
v v ' A au h t dr 

We now develope />, and X, by the powers of 0; so that (S°) and (P°) will 
have only p°, and X", which are elliptic functions of t only, with the constant 
elements a n e n &c. 
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And by substitution, 

s=(S)+0 ___ i( ^^ + ^ dT dr + ^ d ^\i dr j 

+ T"^~ + I( F ^ ) aV~2 {n ~&T + T / dr <l> dr i + h J dr* + d^ 
To this we must add, 

^)-V ) -h"av + T,~Jr~ h, dr 2 *' 



5. In concluding this paper, I will give a brief sketch of a transformation of 
the differential equations which determine the place of a disturbed planet, 
which I have never seen noticed, but which I think is deserving of some con- 
sideration. 

Make Q = — - + R ; the equations with reference to a fixed plane are, 
**.dQ d*y dQ_ **.dQ 

dt i + dx~ Vl di? + dy- U ' d? + ~fc~ (1) 

Suppose now a plane always to make the constant angle i with the fixed 
plane, i being the mean inclination of the orbit to this plane. And let this 
plane slide on the fixed plane in such a manner that their intersection may have 
a uniform motion equal to that of the node, or so that the intersection of the 
two planes may be always the mean place of the node. And let 6 — ant be the 
longitude of this intersection, being constant. 

Let now v t and y, be co-ordinates making the angle 6 - ant with x and y, 
x, lying in the intersection of the sliding and fixed planes. We have 
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x = x t cos (0 — ant) — y, sin (0 — ant), 
y = x, sin (0 — ant) + y, cos (0 — ant). 

We form the two equations, 

_ sm(e _^)^ + ^J + cos (6-^)^ + ^ = 0. 



But 



Therefore, 



. dQ . , n s dQ dQ 

cos (0 — ant) -r- + sin (0 — an£) -=— = -s— , 

fn ,x <*Q ,„ JX dQ dQ 

-sm(e-ant) —+cos(e-ant) j-^-- 

. x d?x . , n .d 2 y dQ - 

(e-ant)- m + sm(B-ant)-£+ : g = 0. 



cos 



df 



df dx, 



- sin (0 - ant) -^ + cos (0 - ant) -£ + ^- = 0. 
Substituting for -^ and -^ , their values in x, and y t , we have 

Or dropping the distinctive marks of x, and y n as tending to produce con- 
fusion, and as no longer necessary : 



d 2 x dQ , _ <fy , . . 



(2) 



We now take the new co-ordinates y t and z t , y, lying in the sliding plane, 
and z t being perpendicular to it ; and, consequently, we have 
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y = y, cos i — z t sin i, z = y, sin i + z, cos i. 

Proceeding by exactly the same steps as before, we find 

d?y, dQ a . dx , , . . . . . x A 

-^ + -j 2an cos % -j — arn 3 cos i (y, cos i — z, sin i) = 0, 

d*V rfQ n ■ -dx „ , . . , . . ., rt 

t^ + j— + 2an sin z -T- + aV sm i (y, cos i — z. sin i) = 0. 
ar dz, dt v,7/ ' 

Or again, dropping the distinctive marks of y t and z, 

d 2 y dQ .dx . . . . • -\ a T 

~ + -; 2an cos « -r — oV cos « ( y cos i — z sin z) = 0. 

ar ay oft ^ ' 

a**z rfQ . a ■ .dx , , . . , . . ., _ 

-7 j + -j- + 2cm sin i -r- + aV sin i (y cos i — z sin i) = 0. 

Again, x, and y t lying in the sliding plane, let x f be the angle tt + Gn£ in 
advance of the line of intersection of the two planes. Then £ may be so de- 
termined, that x t shall lie in the mean place of the apse when projected on the 
sliding plane, or in any other position we please. And as before, 

x — x t cos (it + tnt) — y t sin (tt + €nt), 
y=x, sin {it + €nt) + y t cos (w + Gnt). 

Making for a moment, in order to abridge, 
A =2an -j- — a 2 n 2 x, B = —2ancosi-j a?n 2 cos i (ycosi — zsini) ; 



(3) 



we shall have from the first of (2) and the first of (3), by exactly the same pro- 
cess as in the two former transformations, 

d 2 x cPv dQ 

cos (tt -f Gnt) -j-} + sin (ir + Gnt) -rf + ~ + A cos (-*+ €nt) +B sin (tt + Gnt) = 0. 

d 2 x d 2 v dQ 

- sin (tt + $nt) -j^ + cos (tt + Gnt) -^+-~-A sin (ir + Gnt) +B cos (tt + Snt) = 0. 

Or, putting for -^ and-p their values, 
vol. xxn. i 
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^2 +^r +i '/- u ' 



dfi dx, 

And again dropping the distinctive marks, 

d?x dQ T _ 
^ + ^- + X = °' 

d* 2 + % + 



(4) 



i = -^^ 2ancos 2 g-2ew+2onsin 2 scos2('7r+en#) l+2on -^ sin 2 ^ sin 2 (»+Cn<) - 
#[ o 2 » 2 +e 2 n 2 -2aew 2 cos 2 |-^o 2 w 2 sin 2 i+f^aVsin 2 i-2oen 2 sin 2 | Jcos 2(*+ent)~l+ 
y ( £a 2 n 2 sin 2 i — 2agw 2 sin 2 ~ jsin 2 (w + €nt) + a 2 n 2 sin i cos iz sin (w + Gnt ). 

iV r =-7-j 2£n— 2an cos 2 ^+2awsin 2 ^cos2 (•*+&!*) >— 2an-^sin 2 ^sin2(ir+&tf) 
y[2aen 2 cos 2 |-aV-e 2 n 2 +^a 2 n 2 sin 2 i+^ 2 n 2 sin 2 i-2aen 2 sin 2 |jcos2(-7r+en?)| 

a; (£a 2 n 2 sin 2 i - 2a€w 2 sin 2 1 ) sin 2 (w + grtf) + a 2 » 2 sin i cos u cos (w + e«<). 

The terms in L and -Sf containing z are of the order of the third power of 
the disturbing force multiplied by sin i, and will not be wanted in these cases 
in which i is large, and, I think, cannot be wanted when it is small. Most of 
the other terms are of the second order ; and as i is constant, they are all of a 
very simple character. 

We must add from the last of (3) 



+ 
+ 



dx 
N — 2an sin i -v- + «V sun (ycosi — z sinz). 



(5) 
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The term in N containing z is of the order of the third power of the dis- 
turbing force multiplied by sin 2 i, and cannot, I conceive, in any case be 
wanted. 

"We may transform the co-ordinates x', y 1 , and z\ of the disturbing body 
in like manner to a sliding plane, and to a similar situation on that plane. 
Thus, marking all the corresponding quantities relative to this body with an 
accent, we have 

x' - x\ cos (0' - a'rit) - y[ sin (0' - a'n't), 
y' = x[ sin (0' - a'n't) + y\ cos (^ - a'n't), J = z\ . 
x' l = x' i , y[= y' t cos i' — z' % sin i' , z[ = y' t sin i' — z' % cos i'. 
x' 2 = x' 3 cos (•*' + G'n't) - 2/3 sin (w + g'ra'*) , 
y 2 = a? 3 sin (w' + S'n't) + j/ 3 cos (w 7 + S'ra'J), 4 = z» , 

But we may transform x' t j/, and z' to the plane of the orbit of the disturbing 
body, if we please. 

To find the latitude, marking the letters in the second member with the 
number of the transformation to which they belong, we have 

z = y 2 sin i + z 2 cos i = x 3 sin i sin (v + Snt) + y 3 sin i cos (v + Snt) + z 3 cos i. 

But a being the sine of the true latitude, s the sine of the latitude relative 
to the sliding plane, and/ the anomaly, or the longitude measured from the axis 
of x after the last transformation ; 

z = ra, Xa = r cos/, y 3 = r sin/, and 2 S = r«. 

Therefore, by substitution in the last, and dividing by r, 

a — sin i \ sin (w ■+• Gnf ) cos/M- cos (tt + SnJ) sin/} + s cos i ; 
or, 

o = smi sin (/+ W + £«£) + $ cos i. (6) 

Thus the latitude is very easily found. 

Still marking the letters according to the transformation to which they 
belong, 

3/1 = y t cos i — z 2 sin i = x 3 cos i sin (w + Era*) + j/ s cos i cos (tt + era?) - z 3 sin *. 
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Let v be the longitude on the fixed plane, measured from its intersection 

with the sliding plane, 

y t = r sin v, 

the rest as before ; then, by substitution, 

sin v = QOsi j cos/ sin (w + Gnt) + sin/ cos (ir+£nt)\ —ssini. 

Or, ' 

sin v = cos i sin (/+ •*• + Gnt) — s sin i. ( 7 ) 

This will give the longitude when i is large ; and when it is small, we can 
without difficulty find v itself without finding its sine. 

I must, in conclusion, repeat what I have before said, that I think this trans- 
formation worthy of consideration, 

gcnthwaite hall, 

near Babnsley, Yorkshire, 

Nov. 13, 1848. 



